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Abstract 



. , We propose an efficient reduced-order technique for electronic structure cal- 

^ ■ culations of semiconductor nanostructures, suited for inclusion in full-band 

O _ 

quantum transport simulators. The model is based on the linear combination 

of bulk bands obtained by the empirical pseudopotential method, combined 
> 
0> '_ with the use of problem-matched basis functions numerically generated from 

the singular value decomposition. The efficiency and accuracy of the pro- 

■^ ! posed approach are demonstrated in the case of the dispersion relation of 

O 

cn I hole subbands in an unstrained GaN layer. 
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1. Introduction 

A rigorous atomistic description of the electronic structure in confined 
systems is of crucial importance to understand quantum transport phenom- 
ena in technologically relevant nanostructures [l|. Traditionally, theoretical 



studies of nanostructures have relied on the multiband kp method [2| in the 
framework of envelope- function approximation (EFA) 3|. In this approach 
the electronic states of the nanostructure are expanded in zone-center Bloch 
wavef unctions of the underlying bulk crystal, and the expansion coefficients 
(envelope functions) are assumed to be slowly varying spatial functions. De- 
spite the numerous approximations involved, envelope function approaches 
based on kp models have been widely adopted, mainly due to a fair com- 
promise between simplicity and reliability. However, there is evidence that 
a full-zone description is critical for an accurate quantitative modeling of 
nanostructures. As an example, quantum-mechanical mixing between the 
zone center F and the zone edge X states {T-X coupling) due to the crossover 
from direct to indirect band gap in some cubic III-V systems has a relevant 
effect on electrical and optical properties ^. In general, full- zone electronic 
structure models are attractive because they allow for an atomistic descrip- 
tion of the band structure and a unified treatment of bound and unbound 
states. The latter feature is important in the description of scattering mech- 
anisms between continuum and bound states p, la, [Zj, which play a crucial 
role e.g. in the dynamic properties of semiconductor lasers and in the quan- 
tum efficiency of light-emitting diodes (LEDs). As a topical example, the 
efficiency droop observed in GaN-based LEDs |8| has been attributed by 
some researchers to Auger recombination processes which, due to the large 



energy gap, promote carriers to high energy states in different bands above 



10|. A correct microscopic de- 



the barrier, thus contributing to leakage 
scription of such interband processes would require a full-zone approach. An 
atomistic description of the band structure is also important for the correct 
inclusion of strain and of the strong piezoelectric charges which are predicted 
in GaN-based heterostructures ll| . The conventional EFA approach is obliv- 



ious of the atomic details of heterointerfaces and implies a slow variation of 
the envelope function. While interface-related effects lying outside the scope 
of conventional EFA can be accounted for, in principle, by exact envelope 
function theories such as those developed by Burt 12(1 and Foreman [l3|, the 



description of realistic strain configurations in partially relaxed structures 
and the treatment of defects call for an atomistic approach. 



An alternative to EFA, a 
orbitals (LCAO) method_|14 



ong the lines of the linear combination of atomic 



bulk bands (LCBB) 



17 



15|,ll6|, is obtained by the linear combination of 



181 1 , which avoids the decomposition of the wavefunc- 
tion into envelope functions by expanding the states of the quantum structure 
in terms of the full- zone Bloch eigenstates of the constituent bulk crystals. 
LCBB guarantees that the physical symmetry of the system is preserved, 
and allows for an atomistic description of surfaces, interfaces, and strain 
[l7|. Compared with exact diagonalization techniques, LCBB enables to se- 
lect the physically important bands and wavevector points [l7|. As a result, 
the number of basis functions can be reduced significantly compared to the 
plane-wave basis. Although LCBB has been applied to large scale electronic 



structure calculations 



J], the method is still too computationally intensive 



to be included in carrier transport device simulation codes. This applies in 



particular to nitride-based nanostructures, where charge rearrangement in- 
duced by the presence of externally applied or internally induced fields has 



to be considered for a realistic device description 19| . Ideally, one should try 
to combine a complete quantum-mechanical description with a full-band ap- 
proach, avoiding the computationally prohibitive load of atomistic methods 
and the inherent limitations of the EFA. 

2. From LCBB to LCBB-SVD 

With a view of the above remarks, we present a simple acceleration strat- 
egy, LCBB-SVD, based on the singular value decomposition (SVD). This 
strategy was inspired by a numerical technique proposed by some of the au- 
thors to generate reduced sets of problem-matched basis functions in guided- 



wave finite-element analysis 
tronic structure calculations 



20 



21| and already demonstrated for bulk elec- 
22|. Although this reduced-order technique can 
be applied to any type of nanostructure (2D, ID, or OD), here we restrict 
our attention to 2D systems with a confining potential U{r, z) = U{z) which 
is constant in the r plane normal to the confining direction z. Before pre- 
senting the details of the reduced-order model, we briefiy summarize LCBB. 



Following the notation in 18|, we write wavevectors and reciprocal lattice 



vectors as K = {k,kz), G = {g,gz), respectively. We assume that the 
nanostructure is composed of a single material and we expand the nanos- 



^This approximation is well justified for semiconductor-insulator heterojunctions where 
the band discontinuity is very large, or for the opposite extreme of weakly confining quan- 
tum wells (QWs). For vertical transport across Ill-nitride multiple QWs, where this 
simplification could be questionable, the proposed method can be extended to take into 



tructure eigenfunction ip a.s a. linear combination of bulk states $n'fc'fci of the 
material considered 

"ip = /_^ An'k'k'^^n'k'k'^ (l) 

n'k'k'^ 

where ^n'k'k'^ = Un'k'k'^{i^, z)e^'^ ''^e^'^^^ and Un'k'k'^ is the periodic compo- 
nent, which can be expanded in the plane wave basis set with coefficients 

Bn'k'k-X9i9z) 

Un'k'k'S^, ^) = ^ffYl Bn'k'k'Sa. gz)e^^''e^^'^\ (2) 

The unknown wavefunction ip must satisfy the Schrodinger equation [H + 
U{r,z)]ip = eip, where H is the Hamiltonian corresponding to the kinetic 
energy operator and the periodic crystalline potential. By projecting the 
Schrodinger equation on the generic state ^nkk, , a separate eigenvalue equa- 
tion for each wavevector k in the unconstrained plane can be obtained with 
an appropriate choice of the expansion volume in the reciprocal lattice space 

X E E Ut{K -k. + G.)StCkUGMn'kk, = 

n',k'^ G^ 

e{k)Ankk, 
where L = 27i/Akz is the length of the nanostructure in the z direction, Ak^ 
is the spacing between k^ values, Uxi^z) denotes the ID Fourier transform of 
the confining potential U{z), Gz is a reciprocal lattice vector along kz-, and 
^kkK^z) = {un,K+G:,\un',K') are overlap integrals between periodic com- 



ponents. The wavevector {k,kz) belongs to the first 2D Brillouin zone 18| . 



account the full details of realistic hetero junctions. 



thus \kz\ < Gzml'^-, where G-^m is the magnitude of the smallest 



reciprocal 



lattice vector Gz. The empirical pseudopotential method (EPM) 23| is used 
to determine both the full-band dispersion Ep'^{k, kz) and the Fourier com- 
ponents BnkkXdidz) of the periodic functions of the underlying 3D crystal. 
The expansion coefficients BnkkXOi 9z) are then used to compute the overlap 
integrals S^^^fJy{Gz)- The sum over Gz in (|3]) can be safely truncated to 



include only the terms Gz = (0, ±.Gzni) in addition to Gz = (0, 0) [18 1. 

LCBB leads to an eigenvalue problem of rank UhUk^, where n^^ is the 
number of wavevectors along kz and Ub is the number of bulk bands included 
in the calculation. The number of points n^^ in ([1]) needed to achieve a given 



accuracy is structure-dependent 19|. Moreover, the presence of a confin- 
ing potential implies a lower bound on the number of points for the Fourier 
representation of the potential itself. In general, a few hundred points are 
necessary for an accurate description of the energy dispersion of typical 2D 
systems. In practical cases, the computation time is dominated by the calcu- 
lation of the matrix itself (which requires the calculation of the Bloch waves 
^nkkj rather than by its diagonalization. We will show that the information 
necessary to build the eigenvalue matrix (|3]) can be efficiently extracted from 
a few numerically generated problem-matched basis functions. To this end, 
it is convenient to write the Hamiltonian H in the form used in complex band 



structure calculations [3[, which explicitly displays the kz dependence of H 
(non-local terms and spin-orbit corrections are not included for conceptual 
simplicity) 

HcMk, kz) = Hl^a'kl + i^i,G'^' + ^g,g' (3) 



with 

Hh,G' = ^^G,G' (4) 

Hh,G' = —9^^G,G' (5) 

H^G,G' = |;^(^' + 2^ ■ » + 9')Sg,g' + ViG - G') (6) 

Rather than repeatedly solving the bulk problem for each (fc, k^) point, we 
diagonalize if in a few selected Up points (the expansion points) evenly spaced 
along kz between and Gz„j2. The eigenvectors computed at the expansion 
points are then arranged columnwise in a matrix X with dimensions uq x 
UbUp, where uq is the number of plane waves used in the EPM calculations. 



By applying the economy-size SVD 2J| X = UTiV^ , we obtain a. ug x UbUp 



unitary matrix of left singular vectors U, a n^np x UbUp diagonal matrix S 
with positive elements (the singular values), and a UhUp x ni,np unitary matrix 
of right singular vectors V. The significance of each singular vector in the 
description of the bands considered is measured by the amplitude of the 
corresponding singular value [2J]. Since singular values typically range over 
several orders of magnitude, just a few of them may be needed to obtain an 
accurate description of the band structure. Having selected a suitable lower 
bound for the singular values, the bulk problem can be efficiently solved for 
arbitrary points in the range < A;^ < Czmf^ by diagonalizing the reduced- 
order Hamiltonian of rank nc 

HcG'ik, kz) = Hh,G'kl + Hh,G'k. + H'g,g' (7) 

where Hq q, = WH^ q,U with a = 0, 1, 2 and the columns of U are the ha 
columns of U corresponding to singular values larger than the lower bound. 



The eigenvectors used to build matrix X should be carefully selected in order 
to prevent spurious bands in the desired energy range of the dispersion rela- 
tion. As a general rule, ghost solutions may appear if bands lying close to the 
energy range of interest are not included in X. Having computed the Bloch 
states in the interval < A;^ < Gzm/'^ with the reduced-order model, the 
additional information needed in the secular equation ([3]) can be obtained by 
symmetry considerations. The wavefunction coefficients associated to (fc, k^) 
and (fc, —kz) are related (to within a phase) by 3| 

B^-^{G) = B*(G) (8) 



with K = (fc, —kz) and G = {g, —Qz)- A phase factor may app 

A App. A], 



y depending 



26J. In the 



of the specific choice of the primitive vectors 
reduced-order representation, Eq. (jHj) can be written as B^^ = Tf^^B^j^, 
where Tq q = U'^Tq qU* and Tq q is the matrix that incorporates the swap 
sequence of the G vectors with the appropriate phase factors. Out-of-zone 
states can be constructed using the periodicity condition (strictly valid if the 
basis set is not truncated) 

Bn,K^GAG)=B^^K{G+Gz) (9) 

leading to Bn,K+G, = Tcg'+g^B^^k, with T^g'+g, = U^Tg,g'+gM and 
Tg,g'+g^ = S{G, G' + Gz). Once all the necessary bulk eigensolutions have 
been computed, the overlap integrals Sj^jf'jjf^, {Gz) can be efficiently computed 
by evaluating scalar products in the reduced-order subspace spanned by tJ . 
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Figure 1: Valence subband structure of an unstrained 25 A thick GaN layer with a confining 
potential of 0.1 eV along [001], computed with LCBB (solid lines) and kp EFA (dashed 
lines). 
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Figure 2: Electronic structure of wurtzite GaN computed with EPM (solid lines) and its 
6 X 6 fe-p fit (dotted lines). 
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3. Application example 

Active regions of blue- and green-emitting optoelectronic devices usually 
consist of multiple lll-nitride QWs. As a numerical example, we calculated 
the in-plane hole subband structure of an unstrained 25 A thick GaN layer 
with a confining potential of 0.1 eV along [001]. The LCBB-SVD model was 
obtained by computing the upper six valence bands at rip = 5 expansion 
points, evenly spaced between and G'^.m/2 along k^. The EPM parameters 



were taken from 



27 



28] . with an energy cutoff corresponding to nc = 197 
plane waves. The application of SVD decomposition yields ni,np = 30 singu- 
lar vectors. The corresponding singular values span five orders of magnitude, 
which confirms that the sampling rate in momentum space is sufficient to rep- 
resent the selected valence bands. In order to eliminate the redundancy from 
the basis functions set, the no = 18 singular vectors corresponding to sin- 
gular values having magnitude larger than 1/100 of the dominant one were 
included in U. The reduced-order model was solved in Uk^ =201 points along 
the kz direction between —Czmf^ and Czmf^- The selected problem-matched 
basis functions alone are sufficient to reproduce the LCBB dispersion diagram 
of the structure with excellent accuracy, and no difference can be appreciated 
in Fig.[T] (solid lines) between LCBB and LCBB-SVD results. A comparison 
between Fortran implementations of the standard LCBB and the reduced- 
order technique confirms that the latter is about 10 times faster than the 
former. 

It is interesting to compare the LCBB results with conventional multiband 



kp EFA models. Common kp implementations 
symmetrized version of the wurtzite Hamiltonian 



11 



br GaN are based on the 



29 



30| . which is numeri- 



cally unstable in structures having significant discontinuities of tlie material 
parameters at interfaces 3l|, |32|. Another important issue in EFA is the cor- 
rect ordering of the differential operators. Reliable and spurious-solution-free 
subband structures are calculated here with a finite element discretization in 



33 



13j , ensur- 



34 



35|. 



real space that includes Burt-Foreman operator ordering [12 
ing the ellipticity of the equations in the framework of standard EFA 
Fig. [2] compares the electronic structure of unstrained bulk GaN computed 
with EPM (solid lines) and a 6 x 6 fc-p fitting near F (dashed lines) obtained 
by least-squares optimization. From Fig. [1] (dashed lines), it can be seen that 
Q X 6 k-p models are able to approximate the bound states of the nanos- 
tructure in the limits of EFA. However, simplified approaches, where bound 
levels are treated as 2D states within EFA while the continuum portion of 
the spectrum is described through bulk states, are intrinsically unable to 
provide a microscopic description of capture (continuum-to-bound) and es- 
cape (bound-to-continuum) processes, since initial and final states are not 



treated on equal footing 



Investigations of capture processes in separate 



confinement heterostructures (SCH) based on kp EFA and Fermi golden 



rule have also 



the carriers 37 , 



ed to unphysical results due to the finite coherence length of 
Moreover, a zone-center description of the electronic 



38 



structure is not suitable to describe scattering processes that involve states 
far from F, a notable example being Auger recombination in wide band gap 



semiconductors 



40 



1C|. Although this problem can be alleviated by full- 
4l|, they usually require multiple expansion points 



zone kp approaches ; 

to cover the entire Brillouin zone with a computational cost comparable to 

EPM and possible added complexity related to interpolation issues J42|, |43 |. 



12 



The present LCCB-SVD model allows for a unified full-zone treatment of 
extended (bulk-like) and localized (QW-like) states with a small overhead 
with respect to the EPM problem for bulk semiconductors. LCCB-SVD is 
therefore suited for the evaluation of scattering rates between extended and 
localized states, a crucial ingredient to investigate vertical carrier transport 
across heterostructures. 

Acknowledgments 

The authors would like to thank Dr. Michele Penna for useful discus- 
sions. This work was supported in part by the U.S. Army Research Lab- 
oratory through the Collaborative Research Alliance (CRA) for MultiScale 
multidisciplinary Modeling of Electronic materials (MSME). 

References 

[1] A. Di Carlo, Microscopic theory of nanostructured semiconductor de- 
vices: beyond the envelope-function approximation. Semiconductor Sci. 



Tech. 18 (1) (2003) R1-R31. |doi: 10. 1088/0268-1242/18/1/201 



[2] L. C. Lew Yan Voon, M. Willatzen, The k ■ p Method. Electronic Prop- 
erties of Semiconductors, Springer- Verlag, Berlin, 2009. 

[3] D. L. Smith, C. Mailhiot, Theory of semiconductor superlattice 
electronic structure. Rev. Modern Phys. 62 (1) (1990) 173-234. 



doi : 10 ■ 1103/RevModPhys . 62 . 173 



[4] L.-W. Wang, A. Franceschetti, A. Zunger, Million-atom pseu- 
dopotential calculation of 7-X mixing in GaAs/AlAs superlattices 

13 



and quantum dots, Phys. Rev. Lett. 78 (14) (1997) 2819-2822. 



doi : 10 ■ 1103/PhysRevLett . 78 . 2819 



[5] J. A. Brum, G. Bastard, Resonant carrier capture by semiconductor 
quantum wells, Phys. Rev. B 33 (2) (1986) 1420-1423. 

[6] M. Abou-Khalil, M. Goano, A. Champagne, R. Maciejko, Cap- 
ture and escape in quantum wells as scattering events in Monte 
Carlo simulation, IEEE Photon. Technol. Lett. 8 (1) (1996) 19-21. 



doi: 10. 1109/68.475764 



[7] L. F. Register, K. Hess, Simulation of carrier capture in semiconductor 
quantum wells: Bridging the gap from quantum to classical transport, 
Appl. Phys. Lett. 71 (9) (1997) 1222-1224. 

[8] J. Piprek, Efficiency droop in nitride-based light-emitting 
diodes, Phys. Stat. Sol. (a) 207 (10) (2010) 2217-2225. 



doi : 10 . 1002/pssa. 201026149 



[9] K. D. Chik, B. A. Richardson, On the origin of the carrier leakage in 
GalnAsP/InP double heterojunction lasers, J. Appl. Phys. 67 (5) (1990) 
2660-2662. Idoi: 10. 1063/1.3454781 

[10] F. Bertazzi, M. Goano, E. Bellotti, Numerical analysis of indirect 
Auger transitions in InGaN, Appl. Phys. Lett. 101 (1) (2012) 011111. 



doi: 10 .1063/1 .4733353 



[11] F. Bernardini, Spontaneous and piezoelectric polarization: Basic theory 
vs. practical recipes, in: J. Piprek (Ed.), Nitride Semiconductor Devices: 



14 



Principles and Simulation, Wiley- VCH Verlag, Weinheim, 2007, Cli. 3, 
pp. 49-68. 

[12] M. G. Burt, The justification for applying the effective-mass approxima- 
tion to microstructures, J. Phys. Condens. Matter 4 (32) (1992) 6651- 



6690. doi : 10 . 1088/0953-8984/4/32/003 



[13] B. A. Foreman, Elimination of spurious solutions from eight- 
band k ■ p theory, Phys. Rev. B 56 (1997) R12748-R12751. 



doi : 10 ■ 1103/PhysRevB . 56 . R12748 



[14] J. C Slater, G. F. Koster, Simplified LCAO method for the periodic 
potential problem, Phys. Rev. 94 (6) (1954) 1498-1524. 

[15] A. Muiioz, I. Nieves, E. Suarez, BANDAS: a program for teaching band 
theory (LCAO) in solid state physics, Comp. Phys. Comm. 121-122 



(1999) 706. |doi: 10. 1016/30010-4655(06)70110-2 



[16] S. Schulz, D. Mourad, S. Schumacher, G. Czycholl, Tight-binding 
model for the electronic and optical properties of nitride-based 
quantum dots, Phys. Stat. Sol. (b) 248 (8) (2011) 1853-1866. 



doi : 10 ■ 1002/pssb . 201147158 



[17] L.-W. Wang, A. Zunger, Linear combination of bulk bands 
method for large-scale electronic structure calculations on strained 
nanostructures, Phys. Rev. B 59 (24) (1999) 15806-15818. 
doi : 10 ■ 110 3/PhysRevB . 59 . 15806[ 



[18] D. Esseni, P. Palestri, Linear combination of bulk bands 
method for investigating the low- dimensional electron gas in 

15 



nanostructured devices, Phys. Rev. B 72 (16) (2005) 165342. 



doi : 10 ■ 1103/PhysRevB . 72 . 165342 . 



[19] F. Chirico, A. Di Carlo, P. Lugli, Efficient self-consistent pseudopoten- 
tial calculation of nanostructured devices, Phys. Rev. B 64 (4) (2001) 



045314. doi : 10 . 1 103/PhysRevB . 64 . 045314 



[20] F. Bertazzi, O. A. Peverini, M. Goano, G. Ghione, R. Orta, R. Tascone, 
A fast reduced-order model for the full-wave FEM analysis of lossy in- 
homogeneous anisotropic waveguides, IEEE Trans. Microwave Theory 



Tech. MTT-50 (9) (2002) 2108-2114. |doi : 10 . 1109/TMTT . 2002 . 8023231 

[21] F. Bertazzi, G. Ghione, M. Goano, Efficient quasi-TEM frequency- 
dependent analysis of lossy multiconductor lines through a fast reduced- 
order FEM model, IEEE Trans. Microwave Theory Tech. MTT-51 (9) 



(2003) 2029-2035. jdoi : 10 . 1109/TMTT . 2003 . 8158751 

[22] M. Penna, F. Bertazzi, M. Goano, A reduced-order technique for 
the acceleration of electronic structure calculations, in: 14th Interna- 
tional Workshop on Computational Electronics, Pisa, 2010, pp. 65-68. 
doi : 10 ■ 1109/IWCE . 2010 . 56779411 



[23] M. L. Cohen, J. R. Chelikowsky, Electronic Structure and Optical Prop- 
erties of Semiconductors, Springer- Verlag, Berlin, 1988. 

[24] G. H. Golub, C. F. Van Loan, Matrix Computations, 3rd Edition, John 
Hopkins Press, Baltimore, MD, 1996. 

[25] M. Hjelm, Monte Carlo simulations of homogeneous and inhomogeneous 



16 



transport in silicon carbide, Ph.D. thesis, Kungl Tekniska Hogskolan, 
Stockholm (Feb. 2004). 

[26] K. F. Brennan, E. Bellotti, M. Farahmand, H.-E. Nilsson, P. P. Ruden, 
Y. Zhang, Monte Carlo simulation of noncubic symmetry semiconduct- 
ing materials and devices, IEEE Trans. Electron Devices ED-47 (10) 
(2000) 1882-1890. 

[27] M. Goano, E. Bellotti, E. Ghillino, G. Ghione, K. F. Brennan, Band 
structure nonlocal pseudopotential calculation of the Ill-nitride wurtzite 
phase materials system. Part I. Binary compounds GaN, AIN, and InN, 



J. Appl. Phys. 88 (11) (2000) 6467-6475. |doi: 10. 1063/1 . 1309046 



[28] E. Bellotti, F. Bertazzi, M. Goano, Alloy scattering in AlGaN and 
InGaN: A numerical study, J. Appl. Phys. 101 (12) (2007) 123706. 
doi: 10. 1063/1.27483531 



[29] S. L. Chuang, C. S. Chang, k ■ p method for strained wurtzite 
semiconductors, Phys. Rev. B 54 (4) (1996) 2491-2504. 



doi : 10 ■ 1103/PhysRevB . 54 . 2491 



[30] S. L. Chuang, C. S. Chang, A band-structure model of strained 
quantum-well wurtzite semiconductors. Semiconductor Sci. Tech. 12 
(1997) 252-263. 

[31] F. Mireles, S. E. UUoa, Ordered Hamiltonian and match- 
ing conditions for heteroj unctions with wurtzite symmetry: 
GaN/Al^Gai_^N quantum wells, Phys. Rev. B 60 (1999) 13659- 



13667. |doi : 10 ■ 1 103/PhysRevB . 60 . 13659 

17 



[32] F. Mireles, S. E. Ulloa, Strain and crystallographic orientation effects 
on the valence subbands of wurtzite quantum wells, Pliys. Rev. B 62 



(2000) 2562-2572. doi : 10 . 1 103/PhysRevB . 62 . 2562 



[33] B. A. Foreman, Effective- mass Hamiltonian and boundary conditions 
for the valence bands of semiconductor microstructures, Phys. Rev. B 



48 (1993) 4964-4967. |doi : 10 . 1103/PhysRevB . 48 . 4964 



[34] R. G. Veprek, S. Steiger, B. Witzigmann, Rehab le k ■ p band structure 
calculation for nanostructures using finite elements, J. Comp. Electron. 



7 (2008) 521-529. jdoi : 10 . 1007/sl0825-008-0262-y 



[35] R. G. Veprek, S. Steiger, B. Witzigmann, Operator ordering, el- 
lipticity and spurious solutions in k ■ p calculations of Ill-nitride 
nanostructures. Opt. Quantum Electron. 40 (2009) 1169-1174. 
doi : 10 ■ 1007/sll082-008-9259^ 



[36] F. Rossi, Coherent phenomena in semiconductors. Semiconductor Sci. 
Tech. 13 (1998) 147-168. 

[37] M. Mosko, K. Kalna, Carrier capture into a GaAs quantum well with 
a separate confinement region: comment on quantum and classical as- 
pects. Semiconductor Sci. Tech. 14 (1999) 790-796. 

[38] N. A. Zakhleniuk, C. R. Bennett, V. N. Stavrou, M. Babiker, B. K. 
Ridley, Quantum capture of injected electrons in GaN-based laser het- 
erostructures, Phys. Stat. Sol. (a) 176 (1999) 79-83. 

[39] S. A. Levetas, M. J. Godfrey, Calculation of capture of carriers by quan- 
tum wells, Phys. Rev. B 59 (15) (1999) 10202-10207. 



[40] F. Bertazzi, M. Goano, E. Bellotti, A numerical study of Auger re- 
combination in bulk InGaN, Appl. Pliys. Lett. 97 (23) (2010) 231118. 



doi: 10. 1063/1.3525605 



[41] R. Beresford, Full-zone k ■ p method of band structure calculation for 
wurtzite semiconductors, J. Appl. Phys. 95 (11) (2004) 6216-6224. 



doi: 10. 1063/1. 1713043 



[42] C. Persson, C. Ambrosch-Draxl, A full-band FPLAPW + k ■ p-method 
for solving the KohnSham equation, Comp. Phys. Comm. 177 (2007) 



280-287. doi : 10 . 1016/j . cpc . 2007 . 02 . Ill 



[43] A. Marnetto, M. Penna, M. Goano, An accurate dual-expansion-point 
fuU-Brillouin-zone k-p model for wurtzite semiconductors, J. Appl. Phys. 



108 (3) (2010) 033701. |doi : 10. 1063/1 .3459883 



19 



